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ABSTRACT 


The steady state Navier-Stokes equations (NS) with slip 
and temperature jump conditions, with and without ablation 
effects on the surface are solved for a stagnation point 
merged layer flow on the spherical nose of a space vehicle. 
Local similar solutions have been obtained and the resulting 
ordinary differential equations (ODE) have been solved by 
method of Successive Accelerated Replacement (ASR) . 

Particular emphasis is placed to understand the physics 
of the real gases through the concept of effective gamma (y) . 
The ambient fluids in different tunnels have different fluid 
characteristics which have widely different values of y and 
Pr. A detailed study of varying y , Pr, surface temperature 
and viscosity-temperature relations has been carried out in 
order to enable the experimentalist to interpret his data 
appropriately in relation with the available results of ideal 
fluids (y = 1.4) . 

Computations have also been carried out for the merged 
layer with ablation, represented by the oozing out of the 
gases from an impervious surface in the presence and absence 
of surface slip and temperature jump conditions. Slip 
conditions modified due to injection of gases from the surface 
have been used. General interesting features of the flow field 
have been pointed out in the text. 
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SECTION I 


INTRODUCTION 


1.1 General Remarks 

In the present investigation, an attempt has been 
made to understand the basic physics of the flow field 
around a space vehicle under conditions that it encounters 
during the early stages of its reentry into the earth *s 
atmosphere. Due to the rarefaction of the flow, the 
shockwave like structure in the outer shockwave like 
region of the disturbed zone of the space vehicle is thick 
and merges smoothly with the viscous flow underneath. Also, 
due to lack of enough number of collisions among the fluid 
particles in the disturbed zone, the particles reach the 
surface with sufficient momentum and slip along it. Conse- 
quently, the temperature at the base of the viscous layer 
is also different from the surface temperature. This 
entire zone of disturbance (including the shockwave like 
structure) with slip and temperature jump conditions is 
termed as Merged Layer (ML) . The full Navier-Stokes (NS) 
equations with surface slip and temperature jump conditions 
are integrated from the surface to the free stream. Shock- 
wave like structure comes out as part of the integration 


domain. 
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As the vehicle descends through the atmosphere, a 
denser region is encountered by the space vehicle and 
consequently, the shock and viscous layer zones become 
thinner. In particular, the incipient merged layer, 
viscous shock layer, vorticity interaction regime and 
ultimately the boundary layer regime pass over the 
stagnation region of the flow field (Ref. 3). 

Here, we strive to understand the various pheno- 
menon associated with the rarefaction effects of the flow 
field. Particular attention has been paid in understanding 
the effect of variation in the wall temperatures, with and 
without surface slip and temperature jump conditions, on 
the flow field structure and overall characteristics (C H 
and C ) - of the flow field. A very interesting feature is 

F T w 

observed. For a very cold wall case (— — <. 0.15), slip incre- 

o°° 

ases the heat transfer and skin-friction coefficients 

in comparison to the noslip case while for hot wall 
T 

( T w - >L 0.15) slip decreases the heat transfer and skin 
o“ 

friction coefficient. Also, we find that for low Reynolds 
number, adiabatic wall temperature is greater than the 
stagnation temperature, while for higher Reynolds number, 
adiabatic wall temperature and stagnation temperature tend 
to become equal . A possible explanation of this phenomenon 
is given in Section - III on "Results and Discussion". A 
thorough investigation of the various viscosity-temperature 
(li- T) relations on the flow field was carried out. It 
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was found that for prescribed free stream conditions, data 

for the heat transfer coefficient against the free stream 

Reynolds number with slip fell on different curves for 

different y-T relations. The different curves collapse 

into a single curve when the same set of data are plotted 

P ~ U “ r B, 

against the stagnation Reynolds number (Re = — ) . We 

^0“ 

call this curve as universal curve in the sense that it is 
valid for all values of stagnation Reynolds number. 

When the space vehicle enters the earth's atmosphere, 
intense heat is generated in the disturbed zone of the flow 
field and the molecules dissociate and ionise. Thermo- 
dynamic properties which are constant for a perfect fluid 
no longer remain constant. Also, to analyse the disturbed 
zone of the flow field, we have to solve species conservation 
equation in addition to the mass, momentum and energy equations. 
Chemical kinetics come into picture which makes the analysis 
all the more difficult, time consuming and tedious. Here, . 
we have tried to incorporate the real gas effects in a simple 
way by using the concept of effective 'y' (Ref. 4). The 
results predict the behaviour of the flow quantities and 
overall flow characteristics only qualitatively - 

This kind of study is extremely useful for wind tunnel 
data as well. In wind tunnels, different types of gases are 
being used which have different values of y and Pr. Nitrogen 
(diatomic) , Argon or Helium (Monatomic) and air are most 
extensively used in the experiments. In NASA Langley Research 
Centre, Freon (with y “1.28) is being used to simulate real 
gas effects in the wind tunnel. 
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The present study helps to understand the effect 

of these gases (with different values of gamma) on the 

flow characteristics in comparison to air (y = 1.4). 

Here, intensive investigations have been carried out and 

the values of C^, and flow parameters are plotted 

against gamma and Prandtl number. Particular attention 

has been paid to study the effect of gamma on the wall 

pressure and the overall flow characteristics at different 

values of • stagnation Reynolds number. It is observed that 

the results for wall pressure are in agreement with the 

(17) 

experimental results of Potter and Bailey 

The space vehicle surface is made up of composites 
like carbon-carbon, carbon phenolic or carbon epoxy. 

During the reentry of the space vehicle, due to intense 
heat in the disturbed zone of the flow field, a chemical 
degradation of these materials takes place. The surface 
gets charred, and the material melts, vaporises or sublimates 
The pyrolysis gases which ooze out mix with the main flow and 
affect the flow parameters and overall flow characteristics. 
In the process of sublimation the space vehicle surface may 
receed. 

In the present analysis, the effect of receeding 
of the surface is ignored- Here, an attempt has been made 
to study the effect of ablation in a very simple way by 
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simulating the ablation process through the emission of 
gases from the impervious surface. The high heat capacity 
of the pyrolysis gas, surface recession and roughness due 
to charring have been ignored. Under such surface condi- 
tions, surface slip and temperature jump conditions are 
modified due to injection of gases in the Knudsen layer. 
The modified slip conditions at the surface are used in 
the integrations of the governing equations . Extensive 
computations under varied . small emission rates of 
gases have been carried and the results are plotted as 
a function of mass transfer parameter, M (M = . . 

r'oo CO 

1 . 2 Review of Literature 

Different approaches have been used to analyse the 

rarefied hypersonic flow past a blunt body. During the 

reentry of the space vehicle, it has to pass through 

various flow regimes which have been classified by 

(3) 

Probstem et al . The flow past the space vehicle 
cannot be solved by one method for all the regimes . 
Consequently different procedures have been adapted to 
investigate these regimes. 

Van Dyke developed second order boundary 

layer theory to evaluate the viscous hypersonic flow 
past the nose of a blunt body on the basis of asymptotic 
expansions of the NS equations. In the process he 
identified seven second order effects: longitudinal 
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curvature/ transverse curvature, slip, temperature jump, 
entropy gradient, stagnation enthalpy gradient and 
displacement. He defined a parameter e (viscous hypersonic 
similarity parameter) given by 

o 2 

[ <Y- 1 ) K ] 

e = 

/~ Re 

CO 

where is the index in p-T law. Van Dyke^ 2 ^ carried out 
the analysis by expanding the flow variables asymptotically 
in powers of e for the inner (viscous flow) and the outer 
(inviscid flow) regions of the flow field. Matching the 
two expansions gave a set of equations valid for the entire 
region behind the shock. Van Dyke^ 20 ^ calculated the flow 
upto second order term. Kao^"^ developed third order 
boundary' layer approximation. He too used the method of 
matched asymptotic expansions . He termed the boundary 
layer and the shock wave as the "inner" layer (viscosity 
is significant in these regions) and the nearly inviscid 
region between them as the "outer" region. After 
expanding the flow variables in the various powers of e 
in the inner and outer regions/ te collected like powers 
of e which yielded a set of equations to third order. 

These equations were found to be valid in the entire disturbed 
zone. Viscous effects appear in the third order equations 
for inviscid flow region. This procedure is generally applied 
to the incipient merged layer regime. The higher order 
boundary layer method gives accurate results only upto 
slightly rarefied regimes. Also, this method is very time 
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consigning. 

Another way of attacking the problem is through 

( 2 ) 

viscous shock layer approximation. Davis developed 

viscous shock layer (VSL) equations to solve the flow 

past axisymmetric blunt body moving at hypersonic 

speeds. These equations are truncated NS equations 

retaining terms upto second order in the inverse square 

root of the Reynolds number. These equations along with 

surface slip and shock slip conditions can describe the 

entire disturbed zone. Implicit finite difference 

( 2 ) 

method was used by Davis to solve these equations. 

Probstein and Kemp also used viscous layer and 

incipient merged layer approach to solve the flow past 
a sphere and a cylinder. The shock layer region is 
considered to be fully viscous. Rankine-Hugoniot 
relations give the outer boundary conditions for the VSL 
approach. Locally similar solutions were used by 
Probstein et al. to estimate the flow in the stagnation 

region . 

The above approaches , though simple to use, are 
valid in the range from intermediate to very high Reynolds 
number. This formulation fails at sufficiently low Reynolds 
number where the other terms that have been neglected 
in the NS equations become significant. 
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Hypersonic Merged Layer (HML) approximation is 

another method to estimate the hypersonic flow past a 

blunt body. This is the only method which holds promise 

and gives good results even at sufficiently low Reynolds 

number and the surface quantities obtained from this 

approach and kinetic approach are found to be in good 

(12) 

agreement. Levinsky and Yoshihara were the first to 

use this method on the basis of thin layer approximation 
with noslip surface conditions. Their results predicted 
a decrease in impact pressure with decrease in Reynolds 

number which was not in agreement with the experimental 

(17) (12) 

results of Potter and Bailey . Levinsky and Yoshihara 

attributed this anomaly to the neglect of slip boundary 

conditions. But later it was shown by Jain and Admurthy ^ 

that inclusion of slip effects further decreased the impact 

pressure. Jain et al . ^ improved upon the analysis of 

( 12 ) 

Levinsky and Yoshihara using full NS equations with slip 

and temperature jump effects to describe the flow field. 

They developed local similar solutions valid in the stagna- 
tion zone of the flow past the nose of a space vehicle. 

(5) 

Jain and Prabha compared the results of HVSL, 

HML and DSMC (Direct Simulation Monte-Carlo Technique) 
with the experimental results to find the limit of validity 
of HVSL theory. It was found that HML gave good results 
upto Kn <_ 0.6 while HVSL gave accurate results only upto 


Kn < 0.06. 
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• ( 9 ) 

Jain developed two term indented series solution 

of the full NS equations for the merged layer flow on the 

spherical nose of a body incorporating the downstream 

( 8 ) 

effects- Jain and Prabha derived the equations governing 
the first two indented terms of the series expansion with 

( 7 ) 

asymmetry effects. Jain v in his series of lectures on 
reentry aerodynamics tried to generate an understanding of 
the flow under hypersonic rarefied conditions for adiabatic 
and cold wall cases . 

Till now, various numerical techniques have been used 

(12) 

to solve the equations of motion. Levinsky and Yoshihara 
integrated the equations from each end and matched the 
solution at a common point between the shock and the 

boundary layer where the integration procedure was stable. 

. ( 9 ) „ 

Jam used the method of Accelerated Successive Replacement 

(ASR) to integrate the equations numerically. In the present 

investigation, ASR method has been used. 

The problem of ablation involves slip boundary 

conditions at the surface with mass transport effects. 

Whitehead and Davis usec j the analysis of Shidlovsky 

( 18 ) 

and Patterson v to give a set of such equations which 
allow for the mass transport across a porous surface. 

( 14 ) 

Lewis et al. solved the full VSL equations for hyper- 

bloids at zero lift and spherically blunted cones at 

large angles of attack. They modified the boundary 

( 22 ) 

conditions given by Davis for the case of ablation. 
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assuming porosity at the wall to be zero. Moss showed 

the effect of mass injection on flow properties within and 
downstream of the injection region under noslip surface 
conditions . 
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SECTION II 


MATHEMATICAL FORMULATION OF PROBLEM 


2.1 Discussion of Governing Equations 

The steady state Navier-Stokes equations in spherical 

polar coordinate system with axial symmetry for an ideal 

(23) 

compressible fluid are given in Yuan ' and (Ref. 8). These 
are as follows : 


Continuity Equation : 


■ (pv) r + + i (pu)<j + pu Co^O 


= 0 


( 1 ) 


Radial Momentum Equation : 


P r + P (v v r + p v 0 - ^-) 


r 


(r v + 2v + u Q + u Cot 0) ] 


= 2 L v r - ~ T ' -• ■ 


+ t r <r> r + J o + 


r 

2u u cot e, 

— (2v r - - I 


+ ^ {Cot o[r(p) r + — ] - — (v + u Cot Q) } (2) 
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Transverse Momentum Equation: 


/ u u V N 

+ P ( v u r + r u © + — } 

~ r[ r (u © + v) “ 3F v r + 2v + u @ + u Cot 0) ] 

3P 


+ (v[ r( f ) r + + F 1 + 3 


O TJ 

+ Cot © (u - u Cot ©) 

r y 


© 


(3) 


Energy Egua.tion : 


p (V h r + a h Q ) 


= v p r + f p 0 


+ Fi** (y r Vr + r h e + yh r + Vh G 


+ y[( 2v r + — 2 ( U G + v) 2 + --j (v + u Cot Q) 


+ r 


C (f) + ] 


2y 

3r' 


[r v r + Uq + 2v + u Cot ©] 


(4) 


Eguation of .State : 

p = pRT 

and viscosity temperature relation is defined as 

U = y(T) 


(5) 


( 6 ) 


Here r represents the radial distance and © represents the 
vectorial angle; u and v represent the tangential and normal 
velocity components in increasing r and © directions (See 
Fig. A); p , p, h are the density, pressure and specific 
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enthalpy respectively. Suffixes r and 0 denote the derivatives 
of the dependent variables with respect to ' r ' and 'O' 
respectively . 


Boundary Conditions 

At the edge of the merged layer r = r & 

u = U Sin © 

00 

V = -U Cos 0 

00 

p = Poo (7) 

P = Poo 

T = T « 


(a) At r = r B# i.e. at the surface of the body, slip 
velocity and temperature boundary conditions are 


u = 


2-o 


E x <fr>l r - v = 0 

B 


T = T w + 


2-g 2y r 
a Y+l L 


L- 111 

Pr 9r /J r 


hi 


( 8 ) 

(9) 


B 


A w is the mean free path of the gas at the body 
surface and 

X 


-/ 


TTY 

2 


pa 


For adiabatic wall case. 


<£>r 


( 10 ) 


B 


(b) 


For the noslip case at r = r 


u = v = 0 and T = 


B 

T w ° r 


= 0 
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2 . 2 Non-Dimensionalisation Procedure 

The variables have been non-dimensionalised as follows: 


u - 


u_ 

U 


v = 


v_ 

U. 


P u 2 

CO CO 


p_ 

p 


( 11 ) 


- Qoo 


y = 


y 


Ooo 


r = 


• r. 


Re = 


B 


U P 

CO g oo 


QOO 


where suffix 1 00 • denotes the value of variable at the free 
stream location and the suffix , o“' denotes its value at the 
free stream stagnation conditions; r fi is the nose radius. 

The governing equations (1) to ( 6 ) in the non-dimensional 
form become 


Continuity Equation: 

( o V) + - 2 -£* - 

r r 


(pu) 


G + p 5 Cqt_Q = Q 


(12) 


Radial Momentum Equation : 


P (v v _ + ~ v - 3 ) 


-2 

u 


© 


Re 


[ u v_ - ( r v_ + 2v + u + u Cot ©)] 


3r 


+ ib E {5 + i a]> e + k J L 2v r 


u_ + V 

!: [2v r -2—.] 


+ ~ {Cot © [r (—) _ + z (v + U Cot ©) } 

“ JT XT IT 3T IT 


( 13 ) 



h I |m 


Transverse Momentum Equation 


P G + P (v u_ + 3— u + 

r r r 


~ - [y (*u + v) 


Re - Lr “9 


(r V + 2 v + 
3 r r 


r r r 


+ |b [ 5 ^ [r <§> + 

Re r r r 


2 y 


Re r 


_2 Cot O (Uq - u Cot ©) 


Energy Equatio n 


A p ( v T_ + — T ) = v p_ + 3- p 
r r w r r 


+ 


Re 


i— =■ f (y r 
Pr r L v M 

T_) 

r 

+ 

r 

T 

- 9 

r 

- r -2 
y [2 v_ + 

r 

2 

_2 

r 

(5 © 

- 2 
+ v) 

Z) _ + zr 

r r r 

i 2 

] 


1 “r P v_ 

: ^ 

+ 

u © 

+ 2v 


U Q +u Cot 



r 


( 14 ) 


T 0 Cot © 
r 

Cot 9 ) 2 


( 15 ) 


Equation of State 
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where 


A 


C T 

P Qoo 


U 



1 

(Y- 1) M 2 
' 00 


The boundary conditions in the non-dimensionalised form 
as follows: 

At r = r 

e 

u = Sin 0 
v = -Cos 0 



y is the specific heat ratio and is the free stream 
Mach number . 

At r = 1 


_ / t r:y 

U J 2 


M 

2-0 1 co 


Ooo p/- 


a Re v T 
v = 0 

and the temperature jump condition is 


(ft) - 5) 


3r 


t- — 2Y - 


M 


a y+ 1 p r Re*/ T 


J T o„ 


-Oco p/- 


3T . 

3r 


where 


'W 


w 


L Oco 


(18) 


are 


(19) 


( 20 ) 


( 21 ) 
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For noslip case 

u = v = 0 and T = T w or ^—3 = 0 (22) 

8 r 

for adiabatic wall case. 

The physical domain of integration is now transformed 
to computational domain by the following transformation 



(23) 


This transformation changes the domain of integration from 
zero to unity. If r eQ denotes the extent of the merged 
layer at the stagnation point of the body, then 



r = 1 + nn eQ 

n eo appears as a parameter in the governing equation. It 

is an additional unknown in the governing equations and is 

evaluated by a criterion ensuring a smooth merging of the 

ML flow with the ambient flow. In the present investigation 

it is found that the velocity profile converges earlier than 

the temperature profile. We have imposed the 

condition that the temperature gradient at r = ? 

e . 

should be negligible. 
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To transform the governing equations into computational 
domain, we use the relation 


3_ _ 3_ In 

3? ^ 3r 


1 


n 


eo 


3 _ 

3n 


( 25 ) 


Consistent with the form of boundary conditions about 

the axis and the fact that Rankine Hugonite pressure behind 

the shock wave and Newtonian pressure on the surface vary as 
2 

Sin ©, we assume the following form of series expansion for 

(»*) 

the various thermodynamic variables 


u 

On# 

©) 

1! 

d i 
0 

(n) 

Sin 

© 

+ U 

2 (n) 

Sin Q + 

V 

(o. 

©) 

II 

<1 

0 

(n) 

Cos 

© 

+ V- 

2 (r l) 

Cos © 

Sin 2 © + 

p 

(tu 

0) 

0 
l a 

II 

(n) 

+ 

^2 

(n) 

Sin 2 

© + .. 

— 

p 

(r). 

©) 

II 

wl 

0 

(n) 

+ P 2 

2 

(n) Sin © + p 4 

(n) Sin 2 © + . 

T 

(n# 

©) 

II 

HI 

0 

(n) 

+ 

f 2 

(n) 

Sin 2 © + . . . 

» m m 

y 

(0, 

©) 

II 

•d! 

0 

(o) 

+ 

V-2 (o) 

Sin 2 © + . . . 

(26) 


Substituting the various expansions of the flow variables 
viz. , equation (26) in the governing equations (12) to (16), in 
boundary conditions (20) to (22), and collecting terms of like 
powers of ©, we get the following set of ordinary differential 
equations and the boundary conditions for the zeroth order 


terms . 
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Continuity Equation : 


1 + n eo r ' 


n 


( p o v o ) + 2p o + V J = 0 


eo 


o v o' 


(27) 


Radial Momentum Equation : 


n 


eo 


izr- + 


0 V v 
P Q O Q 


4 5o n eo 


n 


) - (- 


eo 


+ 


n u 
eo ^o 


2 V ° 
) 


« n 

1 +rirL eo 

eo 


u 


+ (- 


eo 


eo' 


4 


2n_(l+ nn„) n U 


eo K o 


U 4- V 
O O 


2 (1+n n eo ) 1 +1 n eo 


(28) 


Transverse Moinentum Equation : 

- H — 

2p 2 


u o 

~2 

n 

eo 


Re 

p o 


[ 


1 + nn 


v u 
, o o 

+ p ( 

H o v n 


eo 


u (u + v ) 

+ _o_o O } ] 


eo 1 +ti n 


eo 


l n , 


eo 


Ji£L 


U, 


V, 
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In the present computations, square root law and Sutherland 
law, both have been used to relate viscosity with temperature , 


The boundary conditions are: 


At 0 = 1 


3 0 = 1, V o =-l , P 0 = 1. T o = 1 - lj, P 2 = 0 
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At T] = 0 


Slip velocity is given by 
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v_ = 0 


Temperature jump is given by 
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For adiabatic wall case 



0 


(38) 


At 0 = o, for noslip case 


u = v = 0 


dT 

- c 

T w or an 


= 0 


(39) 

(40) 


In the above set of equations (27) to (33) and boundary 
conditions (34) to (40 ) , a prime over a variable denotes the 
derivative with respect to 0 . The ellipic nature of the 
^ equations comes due to the presence of higher order terms, 
'Wiz. u 2 . To. p 2 # P 2 etc.Jaii( 9) has shown that the effect of 
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this second order term near the stagnation region is small 
and these terms can be neglected in comparison to leading 
terms of the expansion in the vicinity of the stagnation 
point- The locally similar solutions transform the partial 
differential equations to a set of ordinary differential 
equations and boundary conditions. These differential 
equations are valid only for small values of ©, i.e. in a 
small region near the stagnation point . 


2 . 3 Concept of Effective Gamma 

i 

Ideal fluids obey the perfect gas law, p = pRT, while 
the real fluids depart from this law and obey p = z'pRT, were 
z* is the compressibly factor. In real fluids the thermo- 
dynamic quantities no longer remain constant. The real 
gas effects include ionization, dissociation and vibration 
degrees of freedom. 


then 


If the density ratio across the shock is given by e 



( 41 ) 


( 4 ) 

From Rankine Hugoniot relation, as given by Probstein , 
we have 



P /P 


(h + e ) - (h + e ) - 
& ® 00 00 



( 42 ) 


Using hypersonic limiting case, i.e. when approaches 
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infinity (M -*<»), in front of the shock, p 

* CO * oo 

tend to zero 
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p s /p s 
h s + e 
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h , e all 

CO CO 


(43) 


The ratio of specific heats, i.e. gararaa, y has no 

significance in the case of dissociating gas. In that case, 

we define Y , "effective ratio of specific heats", and it 
e 

has a value different frora 1.4 for a real gas. 

(4) 

Frora Rankine-Hugoniot relation , in the limiting 

case, M oo , we have 

00 


= I - 1 

lira y + 1 


(44) 


If z denotes the deqree of freedom which depends upon the 

(4) 

motion of the molecule possesses , then 


£ 


-1 

lira 


z + 1 


(45) 


For a monatomic gas, z = 3 (three translational degrees of 
freedom.), and this gives e m = 1/4. Similarly for a 
diatomic gas (three translational and two rotational degrees 
of freedom), e -^im = 1/6- When the diatomic particle is 
excited to vibration, e lim = 1 / 8 . Thus we observe that 
with introduction of real gas effects, e lim decreases, which 
in turn implies from equation (44) that the value of gamma, y 
decreases. In other words, by decreasing the value of gamma, 

Y we can incorporate real gas effects upto some extent. With 
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real gas effects, at elevated temperatures, the value of 
e^irn may be as l° w as °*0 7 which implies a value of 1.15 
for gamma, Y . 

2.4 Mathematical Formulation for Ablation Problem 

The slip flow boundary conditions without mass 

(19) (18) 

transfer were derived by Shidlovskiy and Patterson 

These boundary conditions assume an impermeable surface 
and zero normal macroscopic velocity at the surface. 

They also assume that the mean free path although small is 
large enough so that there is no interaction between the 
incident and reflected molecules. But there occurs a small 
thin layer between the viscous layer and the body surface 
where the interaction among the molecules takes place and 
the velocity distribution of the .molecules is affected. 

This thin layer is called the Knudsen layer (See Fig. B) . 
vj is the injection velocity-. 

To investigate the flow at low Reynolds number with 
ablation, it is necessary to consider boundary conditions 
with mass transport across the surface. The flow near the 
wall is non isentropic and we consider a non-Maxwellian 
distribution in the Knudsen layer. 

If F denotes the Maxwell’s velocity distribution function, 
then the non-Maxwellian distribution function in the Knudsen 
layer is given by F s , where 
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and 

G is a 


The slip flow boundary conditions given by 

(19) (22) 

Shidlovskiy were modified by Whitehead and Davis 

for mass transport effects by assuming specular and diffused 

(14) 

reflections at the surface. Lewis used the same equations 

for calculating the HVSL flow with ablation. The various 
expressions for the transport of mass, momentum and energy 
across the Knudsen layer were derived by considering that 
the net flux at the outer edge of the Knudsen layer equals 
the incident flux, plus the specularly and diffusely reflected 
fluxes at the wall. 

Let the number of diffusely reflected molecules be 
proportional to © and the number of molecules which have 
undergone specular reflection is proportional to (1 - ©) . 

The number of diffusely reflected molecules are further 
subdivided into those that reflect diffusely from the 
solid portion of the surface (proportional to 1-3) and those 
that pass through the pores of the surface and reflect 
diffusely from the interior of the body (proportional to 3) 


F (1 + G) 


(46) 


F = N ( 2irRT) * e 2 s 


-3/2 


lr,2 


(47) 


component of thermal velocities given by Patterson 
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G = ai H ± 4- i H.Hj + | * 1Jk 


(48) 


where 
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Injection area 

Total surface area 


( 49 ) 


For specular reflection 


F s (U, V, W) = F g (U, -V, W) (50) 


The mass transport equation can now be written as 

% 


00 00 OO 00 —V 00 
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— oo o — oo 

: ■ ■ : .. ' : : I' 

' . \ - jV- : v; . . j 

The term on the L.H.S. represents the net mass transported 

■ ' ■ ' . ■ ' j 

at the outer edge of the Knudsen layer across a unit area 

in a unit time. The first term on the R.H.S. gives the 
mass transported across the Knudsen layer by the molecules 
incident on the surface. The second term on R.H.S. represents 
the mass flow of molecules undergoing specular reflection. 
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The third term gives the mass of molecules which are diffusely 
reflected from the interior of the surface with a Maxwellian 
distribution relative to the injection velocity v^ . The 
fourth term gives the mass of molecules diffusely reflected 
from the solid surface with zero macroscopic velocity. 

In the similar way the equations for momentum and 

energy transport across the surface can be derived. as done 

(2 2 ) 

by Whitehead and Davis . After non-dimensionalising these 
equations in the manner shown in equations (11) and reducing 
them by considering small injection velocity (Ref. 22), we 
have 

Mass Transport 

; , 

N s / ^ f s v* +. 2 -3 ^ 

Y S 

= 2(1 -3) N w / ^ T w (52) 


3 (2 +/ 2 ir v.p ) ] 
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Normal Momentum Transport 

P. = BP t [ 1 + 2/ -(v* - v )]+ (1-3) P [ 1 - 2 ] — v* ] 
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i s 


( 53 ) 
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Tangential Momentum Transport 
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Energy Transport 
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In the present investigations, we envisaged that the surface 
material ablates through the efflux of gases from the surface. 
For ablation, we assume the area of the porous wall is zero 
(i.e.3 = 0) . Using this in the above equations we got the 



■following expressions for slip and temperature jump for 
ablation (Ref. 14) : 
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For v = 0, the above conditions reduce to the expression for 
no ablation. 

In the present investigation, surface slip and temper- 
ature jump conditions given by equations (57) and (58) are 
used in the computation of the governing equations (27) to 
(33) . 

2 . 5 Numerical Method of Integration 

The present set of equations consists of three second 
order equations, viz. equations (28) to (30) . The equation 
of continuity (27) is of first order. In the above set of 
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equations , the first order derivatives in the merged layer 

equations , viz. (27) to (31) are replaced by central 

differences and the end point first order derivatives in 

the boundary conditions (34) to (40) are replaced by three 

point one-sided differencing. This way, a second order 

accuracy is achieved throughout the calculations. If X N 

represents the value of any variable (u, v or T) at the 
th 

N mesh point, then 



X N+1 X N-1 
2A 


(59) 


and 
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X N + 1 " 2X N + X N- 


N. 


(60) 


where A = 


and NDIV denote the grid-size and the 


NDIV 

number of divisions into which the merged layer has been 
divided. Here N = 1 represents the “surface and corresponds 
to p= 0 while N = NDIV + 1 represents the free stream 
location corresponding to i) = 1. The resulting algebric 
equations are solved by the method of Accelerated Successive 
Replacement (Ref. 8) . 

Using the finite differencing the boundary conditions 
(34) to (40) become 
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v 0 (1) 


(62) 
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U Q (NDIV + 1) = 1, V Q (NDIV + 1 ) = -1, f Q (NDIV + 1) = 1 - 


P Q (NDIV+1) = 1, p 2 (NDIV + 1) = 0 


(64) 


For no slip case 


u o (1) = v o (1) = °* T o (1) = T w 

U Q (NDIV+1) = 1 , v Q (NDIV+1) = - 1 , T q (NDIV+1) 
p Q (NDIV + 1) = 1 , p 2 (NDIV +1) =0 


(65) 


1 - 
( 66 ) 


1 
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Accelerated Successive Replacement Method (ASR) 

(13) 

The ASR method, originally developed by Liberstein v , 
has been applied by Jain and Adimurthy ^ to hypersonic 
flows in stagnation region. They found the results to be 
reasonably accurate for a wide variety of parameters. As 
the name implies, the scheme replaces the old values of the 
variables as soon as the new values are available. Also, the 
correction procedure involves a parameter * oj * known as 
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the acceleration parameter* The way this parameter is evalu- 
ated ensures that the correction applied to the variables is 
always small. As such, the iteration procedure does not 

diverge even' under quite stringent conditions. The correction 

*t In t "h 

to a variable at the N mesh point in (k+1) iteration is 
applied in the following manner: 



th 

where x N denotes any variable u, v, T at the N grid point 

and X N (u Q * v Q , n eo ) = 0 is a typical zeroth order 

equation. 'u* is calculated in the following manner to avoid 
divergence of the iterative scheme: 



where e is a small parameter. 

Let 

X N e dX N 

° X N dx N 

If o < l, then we take w = o 
and if & >.1/ then we take w= 1 


( 69 ) 


This procedure ensures that correction applied at each 
iteration is never greater than e times the value of the 
variable at that point in the previous iteration. The 
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correction for each variable starts from N = 2 to N = NDIV. 
In this procedure, boundary conditions are always satisfied 
in each global iterations. 


Solution Procedure 

To start with, an initial profile satisfying the 
boundary conditions and a value of parameter n eQ is pres- 
cribed. Then we proceed to solve the differential equations. 
The corrections are first made to v Q , u Q , T^ for grid points 
2 to NDIV by using equations (28), (29) and (30) respectively. 

These p Q and u 0 are calculated from the equation of state 
and y-T relation respectively , p and p are obtained by 
numerical integration of the first order equations (27) and 
(31) respectively. These equations are integrated from the 
free stream to the surface. At the last but one point from 
the free stream, the values of p Q and p >2 are calculated 
using Trapezoidal rule. The values at the next four subsequent 
points are found by Simpson’s rule and at the rest of the 
points except that on body surface (r)= 0) by Weddle's rule. 
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J Po 


From equation (27) we have 
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It is observed that for v (1) = 0, there occurs a singularity 

o 

and the density at the surface becomes infinite, unless the 
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following condition is used: 

_ 1 __ 

v o — -2n eo U 0 (D for slip case 
= ® for noslip case. 

The above condition is used to improve upon the value of 
v Q at the second grid point. 

Pressure at the surface is calculated from the 
following expression: 

P 0 (1) = P 0 (2) - Ap^l) (71) 


where p Q (1) is calculated from normal momentum equation 
evaluated at the surface, viz. 


P 0 (D 


4 v 0 (D 


3n e0 
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Then knowing the pressure at the wall, the density at the 
wall is evaluated from the equation of state. 

For convergence, we define ERRl and ERR2 as the 

maximum difference in the value of tangential velocity and 

temperature respectively between the current and previous 

tin 

iteration at the N mesh point. When ERRl and ERR 2 are 
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both less than 0.001, the results are considered to be 
converged. The converged results at a particular Reynolds 
number may be used as the initial profile for a higher 
Reynolds number. This reduces the C.P.U. .time significantly. 

However, the following difficulties arose, especially 
at very high Reynolds number: (i) Due to infinitesimally 
thin shock, there are inadequate number of grid points 
in the shock structure zone. For such cases, convergence 
became difficult to attain. (ii) Under such conditions, the 
convergence upto third decimal place was also accepted when 
it was observed that the values of the parameters n , C H 
and (C F /Sin©) have stabilized sufficiently. In a few cases, 
the number of iterati'ons were extended upto 4000 in order to 
ensure that the computations did not diverge. (iii) It was 
observed that the initial profile did have an effect on 
the results. For such cases, the initial profile for computa- 
tion of the results of higher Reynolds number was changed to 
another converged profile at a low Re.(iv) In a few cases, the 
ML thickness- n eo was not permitted to grow. Using this 
technique, we were able to achieve convergence upto fourth 
decimal place in most of the cases and upto third decimal 
place in some cases of very high Reynolds number - 

For ablation also, the above numerical technique of 
ASR method was used. Here, a constant value of n eQ was given 
for both noslip and slip case in order to study the effect of 
ablation on the flow characteristics. Again, for higher 
normal velocity at the surface. We feel that the numerical 
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technique has to be modified if massive blowing is to be 
applied. This may be due to the fact that hypersonic 
rarefied flow is very sensitive to extraneous disturbances 
and the numerical method develops instability when a normal 
velocity is imposed at the surface. 
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SECTION III 
DISCUSSION OF RESULTS 

In the present investigation, an attempt has been made 
to understand the basic nature of the flow around the stagnation 
point of a space vehicle under hypersonic rarefied conditions. 

An effort has been made to imbibe the real gas effects through 

' v 

the concept of effective y (gamma) and to simulate the different 
fluid flow conditions in the wind tunnels. Also, the effect 
of ablation on the flow field is studied in detail. This study 
helps in providing useful information to the designer. 

Calculations have been carried out to study the effect 
of the variation of the surface temperature T^, gamma (y) / 

Prandtl number (Pr) and the viscosity-temperature (y-T) relations 
on the detailed structure and overall characteristics of the 
flow field. In most of the cases, computations have been carried 
out for shuttle reentry conditions at 104.33 Kms and 92.35 Kms 
tfhich respectively represent a rarefied limit of the continuum regime 
and a flow at high Reynolds number. For each prescribed ambient 
rendition given in Table 1, stagnation Reynolds number has been . 
calculated using Sutherland law for viscosity-temperature relation. 

* 

j . i Effect of the Variation .of th e Surface Temperature 

The effect of variation of free stream Mach number (M a ) 
md free stream Reynolds number (Re.) on the flow characteristics 
,„d surface quantities has already been considered by Jain and 
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Adimurthy and Jain and Prabha^^. Their calculations were 
carried out for a very cold wall case and it was observed 
that the heat transfer coefficient, C * with the surface slip 
conditions was more in comparison to the values without the 
surface slip conditions for a wide range of Reynolds number. 

Liu carried out the calculations for a hot wall case and 
found the values of the heat transfer coefficient, 'CL,, 

ri 

without the surface slip conditions to be more than the values 

with surface slip conditions. In the present investigation, 

we have tried to resolve this controversy. In Figs. 1-3, 

we have plotted the variation of the heat transfer coefficient, 

C H , against the wall temperature, T , for the prescribed 

reentry conditions of the shuttle at an altitude of 104.93 Kms 

and 92.35 Kms. From Fig. 1, it was found that for a very cold 

wall case (T <0.15), the surface slip conditions increased 
w 

the stagnation point heat transfer coefficient, C^, while for 
a relatively hot wall case (T 0.15), the surface slip 

conditions decreased the heat transfer coefficient in comparison 
to the' noslip values of C H . 

Stagnation point heat transfer coefficient is given by 



2A y Q 

(1) 

3T 

Pr 

Re 

3? 


For a given ambient condition and Prandtl number, the value 

of C H depends essentially on the product of the viscosity at 

— 3 T 

the surface ( y_(l)) and the temperature gradient (— ) . 

° 3r 

Figs . 2 and 3 respectively show the variation of 
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]i D (l) and ( —) against the wall temperature T . It was 

8r w 

observed that the product of the two quantities, viz. 

— 3t — 

U 0 (l) and (~ for a very cold wall case (T ir <. 0.15) with 
gr^ w 

surface slip conditions happened to be more than the value 

of the product of U (1) and — without surface slip conditions. 

3 r ; 

For a relatively hot wall case (T > 0.15), the value of 

this product with surface slip conditions decreased in 

comparison to the case of no slip conditions. This may be 

the reason for the increase in the stagnation point heat 

transfer coefficient with the surface slip conditions for a 

very cold wall case and decrease in the value of C for a 

H 

hot wall case in comparison to the noslip values of heat 
transfer coefficient. In order to confirm these results, 
we further carried out a similar study for Re = 50 and 
M = 20 (Table 2.2). Fig. 4 shows a similar kind of behaviour 
of the heat transfer coefficient, Cjj, with the variation in 
wall temperature T w . Fig. 5 gives the variation of the 
skin-friction coefficient parameter, C F /SinO, with surface 
temperature T w (See Table 2.1). The skin-friction coefficient 
parameter is also found to behave in a way, similar to the 
variation of the heat transfer coefficient, Cjj, with surface 
temperature T w . It was also observed, that the difference 
in the values of stagnation point heat transfer coefficient, C H , 
and the skin-friction coefficient parameter, C F /SinO, for slip 
and noslip cases becomes small at higher values of the 
Reynolds number in comparison to the difference at lower 
Reynolds number. Figs. 6 and 7 respectively showed the 
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variation of surface slip and the temperature values against 
the wall temperature T . We found that as the wall temper- 
ature increased, the values of the surface slip, u , and 

£z> 

temperature T s , increased for both the altitudes. As 

expected, the values of u o and T_ are higher for a rarer 

s s 

atmosphere at 104.93 Kms than for a denser atmosphere at 
92.35 Kms. We further noted that for the rarer atmosphere 
u £ and (T s - T w ) are not equal to zero even for a very cold 
wall. 

Another interesting feature observed from the above 
analysis (Figs. 1-4) was that the adiabatic wall temperature 
with surface slip conditions was found to be more than the 
free stream stagnation temperature for a rarefied flow 
(104.93 Kms altitude), while for a denser flow (92.35 Kms 
altitude) , the surface temperature tends to become equal 
to the free stream stagnation temperature (See Tables 2.1 
and -2.2) . In the rarer region, due to inadequate number of 
collisions, the molecules strike the surface at an angle 
of incidence (© > 0) and as such carry more translational 
energy E ± to the wall, while the emitted molecules carry 
energy E w corresponding to the wall temperature. The 
energy transferred to the wall at any instant is 
E = E i - E w , which raises the equilibrium temperature. In 
denser region, the molecules suffer enough collisions and 
lose their entire normal momentum by the time they reach 
the wall. Thus, the molecules come at an angle of 
incidence 0=0. As such, they carry a relatively lesser 
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amount of translational energy to the wall, viz., E = E ± -E w 
will be less in the denser medium than in the rarer medium. 
This partially explains the phenomenon observed above. 

3.2 Effect of Different Viscosity -temperature Relations 

From kinetic theory considerations, Bird^ has 

found that the viscosity-temperature variation is of the 

to 

nature of y ’ « T , where the parameter to depends upon the 
model of the gas and intermolecular forces . He further 

(9) 

found, that to = 0.5 for a hard sphere model. Jain 

carried out extensive studies for different to = 0.5, 0.75 

and 1 using the present code and found that for to = 0.5, 

the results of the maximum temperature in the ML agreed 

( 21 ) 

with those of DSMC results of Vogenitz and Takata - In 

the present investigation, the results with to = 0.5 have 

been compared with Sutherland law of y-T relation. The 

(9) 

present results thus supplement the results of Jain 

Figure 8 shows the tangential velocity (u) and 
normal velocity (-v) profiles with the square-root law 
and Sutherland law of y-T relation at y =1.4, Pr = 0.72 
and at an altitude of 104.93 Km. Considerable difference 
was found in the velocity profiles for Sutherland law 
and square— root law of viscosity. Fig. 9 gives the 
temperature profiles obtained from the two viscosity laws 
with and without surface slip and temperature jump conditions. 
In comparison to the square— root law, the results with 
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Sutherland law gave (i) an increase in merged layer 
thickness, (ii) decrease in maximum temperature and the 
shift in maximum temperature towards the outer edge of 
the merged layer and (iii) a diffused shock. These differ- 
ences in the flow characteristics with two y-T relations 
diminished significantly as the ambient atmosphere became 
dense at 92.35 Kms (Fig. 10). From Table 5 (See Section VI), 
we found that the values of the heat transfer coefficient, 

C H , with Sutherland viscosity law are more than the values 
of with square-root law of viscosity by about 16% 
at altitude 104.93 Kms, 26% at 99.49 Kms. and 35% at 92.35 
Kms. Similarly, the skin-friction coefficient parameter, 
C^/SinO, has more value with Sutherland law than with the 
square-root law of viscosity and the percentage difference 
goes on increasing with the increase in Reynolds number. 

Figure 11 shows a smooth curve for the stagnation point 
heat transfer coefficient at different altitudes of the 
shuttle flight conditions for the Sutherland and square-root 
law of viscosity relations. This curve compared well with 
the results obtained from AFE Trajectory ^ 10 ^ and Fay and 
Riddle's empirical formula. We call this curve a universal 
curve in the sense that it is valid for all values of 
stagnation Reynolds number . 
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3-3 Effec t of th e Variation in Gamma and Prandtl Number 
on__the Flow Quantities 

In Figs. 12-28, an exhaustive study has been carried 
out to determine the effect of the variation in Prandtl number 
and gamma on the structure and the overall characteristics of 
the flow field. 

( a ) Effect of Variation of Prandtl Number 

Figure 12 shows the variation of the heat transfer 
coefficient, Cjj, against Prandtl number for different values 
of gamma, y with and without surface slip conditions at an 
altitude of 104.93 Kms of the reentry conditions of the 
space shuttle. It was observed that with decrease in Prandtl 
number, the values of stagnation point heat transfer coefficient 
increased. It was further noticed that as gamma, y decreased, 
the values of the heat transfer coefficient increased for all 
Prandtl numbers. Figs. 11 (b, c, d, e) show the variation of 
the heat transfer coefficient, C H with Prandtl number for 
y = 1.4, 1.3, 1.2 and 1.1 respectively. It was observed that 
the values of the heat transfer coefficient, C^, with surface 
slip conditions were more than the values with noslip conditions 
But with decrease in gamma y, the values of C H with the slip 
and noslip conditions came close to each other mainly 
because °~ — (— °° = l + M ) decreased, which resulted 

T °° T °° 00 T w 

in an increase in the wall temperature ratio T w (T w = - — ) . 

o°° 

With decrease in the Prandtl number, the other changes 
observed in the merged layer were (i) the merged layer thickness 
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increased slightly (Fig. 13), (ii) the maximum temperature 
decreased (Fig. 14), (iii) and there was little effect on 
ths location of the maximum temperature in the merged layer 
under the prescribed conditions (Fig. 15). (Also see 
Tables 3.1, 3.2). 

A similar kind of study was carried out at an altitude 
of 92.35 Kms of the shuttle flight conditions in order to 
confirm the above drawn conclusions. Figs. 16 and 17 respect- 
ively show the variations of the heat transfer coefficient, 

C H' an< ^ *'* ie temperature profiles with different Prandtl numbers 
(See Table 3.3 and 3.4). 

Prandtl number is given by the following relation: 



K 


where k represents the coefficient of thermal conductivity. 
With decrease in Prandtl number, the conductivity increases 
thereby causing more diffusion of heat towards the outer 
edge of the merged layer. As a result, the merged layer 
thickness increases and the maximum temperature decreases. 
Also, larger amount of heat is transferred to the wall, thus 
increasing the heat transfer coefficient, C R , at the wall. 

Figures 18 and 19 respectively show the variation of 
the skin-friction coefficient parameter, C F /SinO, with the 
Prandtl number at altitude 104.93 Kms and 92.35 Kms. As 
expected, the skin-friction coefficient parameter is 
observed to follow the similar kind of variation as shown by 
the heat transfer coefficient, C H . 
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(b) Effect of the Variation in Gamma 

Figures 20-28 show the effect of the variation in 
gamma Y , on the structure of the merged layer and its overall 
characteristics . 

Figure 20 shows the variation of the stagnation point 
heat transfer coefficient , C R with gamma for Pr = 0.6, 0.8, 

1.0 at an altitude of 104.93 Kms, with and without surface slip 
conditions. The figure shows an increase in the heat transfer 
coefficient, C^, with decrease in the value of gamma. Figs. 

20 (b, c, d) show the variation of C„ with gamma for Pr = 0.6, 
0.8 and 1.0 for surface slip and noslip cases. The difference 
between the two curves with and without slip decreased as 
gamma decreased. 

From Figs. 21-23, it was observed that at the altitude 
of 104.93 Kms, the decrease in gamma resulted in (i) decrease 
in the merged layer thickness, (ii) decrease in the maximum 
temperature and (iii) the shift in the location of the maximum 
temperature towards the wall. From Fig. 2 3 we found that the 
effect of the change in gamma was realised more in the shock 
wave like zone of the merged layer, than in the inner viscous 
layer. Figs. 24 and 25 respectively show the tangential 
velocity profiles and the temperature profiles for various 
values of gamma at Pr = 0.8 with surface slip conditions at 
104.93 Kms altitude. In Figs. 26-28, a similar kind of study 
has been done for the reentry space shuttle conditions at 
the altitude 92.35 Kms. The results at high Reynolds number 
(altitude = 92.35 Kms) are similar to the results at low 
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Reynolds number (altitude = 104.93 Kms). However, the 
increase in the value of the heat transfer coefficient, C^, 
with decrease in gamma was observed to be much smaller at 
high Reynolds number than at low Reynolds number. Figs. 29 
and 30 respectively show the variation of the skin-friction 
coefficient parameter, C F /SinO , with gamma at altitudes of 
104.93 Kms and 92.35 Kms. The skin-friction coefficient 
parameter, C F /SinO, is observed to behave in a way similar 
to the behaviour of the heat transfer coefficient ,C H . 

We can physically explain the effect of gamma y , on 
the flow characteristics as follows: 

Gamma, is given by the ratio of the specific heats, 

C 

viz., y = ~ where C and C v are specific heats at 

v p . . 

constant pressure and volume respectively. With decrease in 

gamma, the specific heat at constant volume, C v , increases 

more rapidly than the specific heat at constant pressure, C p . 

For the same energy content (e = C y T) , the maximum temperature 

will decrease with increase in C v . As explained earlier, 

by decreasing gamma, we are to some extent simulating the 

real gas effects which involve ionization and dissociation 

processes. Also, a decrease in y increases density xn 

the merged layer. For the same mass flow rate, the increase 

in density implies a decrease in merged layer thickness due to 

which the heat transfer coefficient, C R , 


increases . 
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3.4 Comparison of the Results from Present Investigation 
with the Available Experimental Data and Theoretical 
Results 

Figures 31-35 show the comparison of the present 
results with the available experimental and theoretical 
results. The overall characteristics of the flow field 
from the present calculations for monatomic and diatomic 
gases have also been compared. 

Figure 31 shows the variation of the heat transfer 
coefficient, C R , with the stagnation Reynolds number for 
different altitudes of the shuttle flight conditions at 
Y= 1.4 and y = 1.66 with surface slip conditions. These 
results were compared with the Direct Simulation Monte- 

( 9 ) 

Carlo Technique (DSMC) results of Moss and Bird and 

with values of the heat transfer coefficient obtained 

from the empirical formula of Fay and Riddle (See Table 4, 

Section VI) . It was observed that the values of C R for 

a monatomic gas (y = 1.66) fell below the values of C R 

for a diatomic gas (y = 1.4). The curve for the heat 

transfer coefficient, C R from the present results at y= 1.4 

(9) 

when compared, with the DSMC results of Moss and Bird 
shows a reasonably good agreement upto an altitude of 
104.93 Kms of shuttle flight conditions. The maximum 
difference between the two results is 8 percent at 104.93 
Kms altitude. The present results are also found to be in 
fairly good agreement with Fay and Riddle values of heat 
transfer coefficient, C H , upto an altitude of 104.93 kms. 
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Figure 32 shows the variation of the skin-friction 
parameter, C^/SinO, for different shuttle flight conditions 
(See Table 1, Section VI) at y = 1.4 and y = 1.66. The 
figure showed that C^/SinO followed a trend similar to 
the heat transfer coefficient, C„ in Fig. 31. 

ri 

In Fig. 33, wall pressure was plotted against stag- 
nation Reynolds number for a monatomic gas (y = 1.66) and 
a diatomic gas ( y = 1.4) at Pr = 0.75 for different 
shuttle flight conditions. It was observed that the wall 
pressure for the monatomic gas had lower values in comparison 
to that of diatomic gas. With decrease in stagnation 
Reynolds number, the wall pressure increased and the increase 
was faster at low Reynolds number. Experimental results 
of Potter and Bailey^ 17 ^ for y = 1.4, T w = T 0& ,M co = 3. 9-6.0 

and for y = 1.66 with T w = T qoo and = 4. 3-5. 7, were also 

(17) 

plotted. Potter and Bailey carried out the experiments 

for a flat-nosed probe and found that the impact pressure for 
monatomic gas dipped to lower values as compared to the 
diatomic gas. The trend shown by the results from present 
investigation compares well with the experimental results 
of Potter and Bailey (17) and Sherman (7) . Fig. 34 shows the 
variation of the wall pressure for different shuttle flight 
conditions at Pr = 0.75 and Pr = 1.0. It was observed that 
the curve for Pr = 1.0 falls below the curve for Pr = 0.75 
and the two curves merged into each other at high Reynolds 
number indicating negligible effect of Prandtl number variation 
on surface pressure in the denser region of atmosphere. 
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3.5 Effect of Ablation 

Figures 35-40 give the affect of ablation on the flow 

variables and overall characteristics of the flow field. 

This study has been done at Re = 50, M ro = 20 and T = 0.2 

with and without surface slip conditions. As pointed out 

in the text, the slip conditions modified for the mass 

transport across the Knudsen layer have been used. Fig. 35 

shows the velocity and temperature profiles with and without 

ablation for surface slip. The surface temperature T g 

decreased with ablation. This is very much evident from 

the relation of T_ from equation (53) . It was noticed that 

s 

the gradients of the flow parameter reduced and merged layer 
thickness increased with ablation. Fig. 36 shows the 
pressure and density profiles. Pressure fluctuations near 
the wall were observed and an increase in the pressure and 
density was noticed in the shockwave like zone of the merged 
layer. The study for the noslip case (Figs. 37 and 38) gave 
a similar kind of variation. Figs. 39 and 40 show a decrease 
in the values of the heat transfer coefficient and skin- 
friction coefficient parameter with increase in the mass 
transfer parameter, H (M = pv) . The present results for 
ablation were found to be consistent with the results of 
Moss (15) who carried out the calculations for the noslip 


conditions . 
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In the present study for ablation, very small 
amount of injection velocity could be achieved (v-j. = 0.002) 
and for larger injection velocities the process becomes 
unstable. The calculations were carried out for a constant 
value of n . A modification in the numerical technique 
for the present code is recommended to apply larger injection 
velocities . 
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SECTION IV 

SUMMARY 


The present study deals with the solution of the NS 
equations to investigate the stagnation point merged layer 
flow past the spherical nose of a space vehicle. 


An attempt has been made to understand the basic 
nature of the flow with variations in the surface temperature, 
viscosity - temperature relations, gamma and Prandtl number. 


Particular emphasis has been placed in simulating the real gas 
effects through the concept of effective gamma, y. Various 
graphs and tables are presented to show the changes in the 
structure and the overall characteristics of the flow field 
under shuttle reentry conditions. 


In the present investigation, we found that slip and 
temperature jump conditions had considerable effect on the 
flow quantities. For a very cold wall case, the surface slip 
condition increased the heat transfer coefficient and the 
skin-friction coefficient parameter from their value with 
noslip conditions while the reverse was true for a relatively 
hot wall. The viscosity-temperature relation was also observed 
to affect the structure and the overall characteristics of 
the flow field. The affect of V-T relation was quite 


significant at low Reynolds number. 


CENTRAL LIBRARY 
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Further, it was observed that a decrease in gamma y, 

resulted in a decrease in the merged layer thickness and 

the maximum temperature in the merged layer while the heat 

transfer coefficient, Cjj increased. A decrease in the 

Prandtl number increased the heat transfer coefficient C , 

and the merged layer thickness but the maximum temperature 

in the merged layer decreased. The comparison of the 

present results with the experimental results of Potter 
(17) 

and Bailey led us to conclude that the surface pressure 

for monatomic gases are lower than the surface pressure for 
diatomic gases. Finally, we found that ablation led to a 
decrease in the heat transfer coefficient, C H , and skin 
friction coefficient parameter, C F /Sin0. Certain modifi- 
cations will be necessary in the numerical technique in 
the present formulation if one wants to apply massive 
blowing. 
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SECTION ut 
T ables 



Y * 1.4, Pr « 0.7 


Altitude p T r D Re M T 
Kras Km/sec gm/crti °K cms °K 

7.47 0 . 115xl0 -9 249.0 136.2 73.148 23.613 387.0 
7.47 0.246xl0 -9 223.0 136.2 171.65 24.95 506.5 
7.50 0.591xl0~ 9 190.0 136.2 474.617 27.14 705.0 
7.50 0.218xl0~ 8 180.0 129.6l750.85 27.88 921.5 
7.53 0 . 637xl0 -8 198.9 132.2 4791.76 26.63 1190.0 


109.75 

104.93 

99.49 

92.35 


85.74 
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2. Values of stagnation point heat transfer coefficient, 

and skin-friction coefficient parameter C^/SinO with surface 
temperature variation 

Table 2.1 


Re' = 

10, M 

oo 

= 20, 

Y = 1.4, Pr = 

= 0.7 





WALL 









TEMPERATURE 








sf 31 

II 

0.05 

0.1 

0.2 0.4 

0.6 

0.8 

0.9 

1.0 

1.1 1.15 

C H 

SOSLIP 

0.6336 

0.643 

0.628 0.5014 

0.368 

0.233 

0.165 

0.097 

0.03 -0.002 

C /SmQ 
F 

0.46 

0.548 

0.608 0.644 

0.669 

0.69 

0.70 

0.711 

0.723 0.728 

i 

C H 

0.644 

0.60 

0.552 0.383 

0.26 

0.144 

0.09 

0.036 

: 1 
1 

-0.013 - | 

SLIP 









C /SinO 
F 

0.5549 

0.558 

0.542 0.547 

0.522 

0.536 

0.535 

0.532 

0.530 - 




I 


Table 2.2 


Re = 50, M to 

It 

>- 

o 

CN 

II 

1.4, Pr 

= 0.7 





WALL TEMPERATURE 

0.05 

0.1 

0.2 0.4 

0.6 

0.8 

1.0 

1.1 

f w 








NOSLIP H 

Cp/Sin© 

0.28 

0.292 

0.275 0.206 

0.139 

0.074 

0.011 

- 0.02 

0.202 

0.228 

0.249 0.262 

0.2*71 

0.279 

0.287 

0.29 

C H 

0.304 

0.295 

0.256 0.198 

0.111 

0.048- 

- 0-012 

- 

SLIP C_/Sin© 

F 

0.23 

0.24 

0.243 0.24380.251 

0.254 

0.256 

— 


3. Values of C H< Cp/Sin©, T 


and n with variation in gamma, 
max eo 


and Prandtl number, Pr 

Table 3 . 1 

Altitude = 104.93 Kms (without surface slip conditions) 


Gamma (- 7 ) 


Pr 


0.6 


H 

Cp/Sin© 

T 

max 


n 


Pr 


0.8 


eo 


H 

Cp/Sin© 


max 


n 


eo 


Pr 


1.0 


'-'H 

Cp/Sin© 

T 

max 


1.4 

1.3 

1.2 

1.1 

0.4557 

0.4991 

0.5606 

0.6752 

0.3352 

0.3675 

0.4141 

0.5178 

0.6011 

0.5814 

0.5579 

0.5528 

0.484 

0.407 

0.326 

0.248 


0.394 

0.44 

0.4915 

0.579 

0.321 

0.356 

0.397 

0.4836 

0.706 

0.684 

0.653 

0.622 

0.457 

0.383 

0.3 

0-215 

0.3584 

0.3954 

0.4539 

0.5508 

0.319 

0.349 

0.397 

0.4652 

0.789 

0.766 

0.734 

0.688 
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Table 3.2 


Altitude = 104. 

93 Kms (with 

surface slip 

conditions) 


GAMMA 

(Y) 

1.4 

1.3 

1.2 

1.1 


g h 

0.704 

0.7047 

0.7102 

0.6814 

Pr=0 .6 

Cj/SinO 

0.546 

0.543 

0.5438 

0.566 


T 

max 

0.648 

0.622 

0.59 

0.572 


n eo 

0.506 

0.423 

0.336 

0.254 


C H 

0.61755 

0.634 

0.64 

0.658 

00 

* 

o 

11 

U 

& 

Cp /SinO 

0.546 

0.549 

0.5488 

0.5573 


T 

max 

0.7579 

0.7298 

0.6913 

0.6477 


n eo 

0.477 

0.397 

0.31 

0.213 


C H 

0.549 

0.568 

0.573 

0.6073 


C F /SinO 

0.55 

0.549 

0.546 

0.557 

Pr=l .0 

T 

max 

0.8401 

0.8121 

0.7728 

0.7183 


rx 

eo 

0.455 

0.379 

0.296 

0.2 
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Table 3 . 3 

Altitude = 92.35 Kms (without surface slip conditions) 


GAMMA (Y) 

see 

1.4 

1.3 

1.2 

1.1 


C H 

0.1956 

0.2174 

0.2474 

0.2919 

Pr=0 .6 

C F /SinO 

0.1235 

0.1322 

0.1441 

0.1642 


T 

max 

0 .9645 

0.963 

0.956 

0.94 


n eo 

0.166 

0.135 

0.102 

0.068 


C H 

0.1655 

0.179 

0.205 

0.2423 

Pr=0 . 8 

C F /SinO 

0.1223 

0.13045 

0.1430 

0.1623 

T 

max 

0.9888 

0.9859 

0.9835 

0.972 


n eo 

0.146 

0.129 

0.097 

0.063 


C H 

C_/ SinO 

0.144 

0.1539 

0.17728 

0.2093 


0.119 

0.12831 

0.1407 

0.16038 

Pr=l .0 

Jb 

T 

1.02 

0.7986 

0.999 

0.9894 


max 

n 

eo 

0.146 

0.126 

0.094 

0.061 
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Table 3.4 


Altitude = 92.35 Kms (with surface slip conditions) 


GAMMA 

(Y) 

1.4 

1.3 

1.2 

1.1 


C H 

0.2315 

0.2466 

0.2639 

0.2932 

Pr=0 .6 

C F /Sin© 

0.1515 

0.1543 

0.16 

0.169 


T 

max 

0.964 

0.9635 

0.958 

0.94 


n eo 

0.1665 

0.1355 

0.1035 

0.0675 


C H 

0.198 

0.2049 

0.2223 

0.2427 


C_ /Sin© 

0.1505 

0.153 

0.158 

0.1685 

Pr=0 . 8 

T 

max 

0.9887 

0.9852 

0.9851 

0.9736 


n eo 

0.1485 

0.1285 

0.0975 

0.0635 


C H 

C F /SinO 

0.176 

0.177 

0.192 

0.2156 

Pr= 1 . 0 

0.1486 

0.1511 

0.1565 

0.168 


T 

max 

1.018 

0.9986 

0.999 

0.987 


n ^ 
eo 

0.1485 

0.1275 

0.0945 

0 .0595 
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4 * Universal C H cur ve ; In this table, the values of the heat 
^" rans ^ er coefficient with surface slip for space shuttle 
reentry conditions (using Sutherland and square— root law of 
viscosity), AFE trajectory and Fay and Riddle's empirical 
formula are compared. Fay and Riddle's formula is given as 
follows : 

rr 1 /A n c 1_ ^ n i 


2 P » U » 


x 0.94) (- 


r t 

“ , 0.5 , qo, 0.1 

n — ) X ( zr — ) 


X (1 + 


i/M x Re 
00 

p 

QA 


(Y-l) m" 


x (1 


where Re is the stagnation Reynolds number 


Y + 1 


•2yM^ - (y-l)' 


l" 1 « [ 


(Y+l) M 0 


-J 

Xri , y T 


(y-l)M +2 


^ (1 + 2 M ~> 


, eo . . Y-l m 2 

and — = 1 + V M a 

00 


Y = 1.4, Pr = 0.72 


Table 4 



SHUTTLE 

VALUES 

FAY & RIDDLE 
VALUES 

AFE TRAJECTORY 

Re 

C H 

C H 

Re 

C H 

4.8132 

0.8479 

0.9774 

10 

0.69 

6.92 

0.7607 

0.8152 

16.6 

0.55 

10.288 

0.6528 

0 .656 

31.62 

0.4 

15.321 

0.5522 

0.5378 

56.23 

0.3 

24.762 

0.43284 

0.41291 

100 

0 .2 

38,99 

0.3419 

0.32906 

158.48 

0 . 145 

139.95 

0.1498 

0.169 

182 

0.13 

258.81 

0 .0998 

0.1222 
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5. Values of the heat transfer coefficient, CL,, and skin- 

ri 

friction coefficient parameter, Cp/Sin© with different 
U-T relations for surface slip conditions 

Table_5 # 

■ Y = 1.4, Pr = 0.72 


Altitude 

Kms 

Square-root law 
of viscosity 

Sutherland law 
of viscosity 

104-93 

C H = 0.5574 

Cp/SinO = 0.45165 

C H = 

C F /Sin0 

0.64824 

= 0.54438 

99.49 

C H = 0.34307 

C F /Sine = 0.259 

C H = 

C„/Sin© 

Jtf 

0.43417 

= 0.3412 

92.35 

C 0 = .0.1498 

H 

Cp/Sin© = 0.10685 

C = 0 

H 

Cp/Sin© 

.20314 

= 0.14727 




FIG. A COORDINATE SYSTEM AND SYMBOLS 
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FIG. B COORDINATE SYSTEM FOR SLIP FLOW BOUNDARY 
CONDITIONS 
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VARIATION OF THE HEAT TRANSFER COEFFICIENT , C H , WITH WALL 
TEMPERATURE FOR SURFACE SLIP AND NOSLIP CONDITIONS 
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\ 3 COMPARISON OF TEMPERATURE GRADIENT WITH WALL TEMPERATURE 
FOR SURFACE SLIP AND NO SLIP CONDITIONS 
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VARIATION OF THE STAGNATION POINT HEAT TRANSFER COEFFIGENT 
C H WITH WALL TEMPERATURE, Tw FOR SURFACE SLIP AND NOSLIP 
CONDITIONS 
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FIG. 5 VARIATION OF SKIN-FRICTION COEFFICIENT PARAMETER 7 C F / sin 6 , 

WITH WALL TEMPERATURE ,Tw, FOR SURFACE SLIP AND NOSLIP 
CONDITIONS 
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FIG. 6 VARIATION OF SUP VELOCITY AT THE SURFACE WITH WALL TEMP- 
ERATURE FOR SHUTTLE CONDITIONS AT 104-93 AND 92-35 Kms. 
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FIG. 7 VARIATION OF THE TEMPERATURE AT THE BASE OF THE VISCOUS LAYER 
WITH WALL TEMPERATURE , Tw FOR SHUTTLE CONDITIONS AT 
104.93 AND 92-35 Kms 




FIG. 8 STAGNATION POINT TANGENTIAL VELOCITY AND NORMAL VELOCITY 
PROFILES FOR DIFFERENT VISCOSITY TEMPERATURE RELATIONS 





FIG. 9 STAGNATION POINT TEMPERATURE PROFILES FOR DIFFERENT 
VISCOSITY TEMPERATURE RELATIONS WITH AND WITHOUT 
SURFACE SLIP CONDITION 
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FIG. 11 COMPARISON OF THE STAGNATION POINT HEAT TRANSFER COEFFICIENT FROM 
THE PRESENT RESULTS FOR REENTRY SPACE SHUTTLE FLIGHT CONDITIONS AND 
AFE FLIGHT CONDITIONS WITH FAY AND RIDDLE VALUES 





Pr.no 


x 

u 


o 


VARIATION OF THE STAGNATION 
WITH PRANDTL NUMBERS AND - 
CONDITIONS 
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POINT HEAT TRANSFER COEFFICIENT ,Ch 
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FIG. 13 VARIATION OF THE MERGED LAYER THICKNESS neo WITH 
PRANDTL NUMBER USING SURFACE SLIP AND NO-SLIP 
CONDITIONS 
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VARIATION OF THE MAXIMUM TEMPERATURE ,Tmax, 
WITH PRANDTL NUMBER AND USING SURFACE SLIP 
AND NOSUP CONDITIONS 
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FIG. 15 STAGNATION POINT TEMPERATURE PROFILES WITH AND WITHOUT SURFACE SUP 
CONDITION FOR DIFFERENT VALUES OF PRANDTL NUMBER 




FIG. 16 VARIATION OF THE STAGNATION POINT HEAT TRANSEER 
COEFFICIENT Ch WITH PRANDTL NUMBER FOR A GIVEN 
VALUE OF GAMMA WITH AND WITHOUT SURFACE SLIP 
CONDITION 
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FIG. 18 
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FIG, 19 


FIG. 18 


FIG. 19 


VARIATION OF STAGNATION POINT SKIN FRICTION COEFFICIENT PARAMETER, 

Cf / Sin 6 ,WITH PRANDTL NUMBER FOR NOSLIP CASE 

VARIATION OF STAGNATION POINT SKIN FRICTION COEFFICIENT PARAMETER, 
Cp/sin 6, AGAINST Pr.no. WITH AND WITHOUT SURFACE SLIP CONDITIONS 
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FIG. 20 VARIATION OF THE STAGNATION POINT HEAT TRANSFER 
COEFEICIENT C H AGAINST GAMMA WITH AND WITHOUT 

SURFACE SUP CONDITIONS 







VARIATION OF THE MERGED LAYER THICKNESS AGAINST GAMMA 
WITH AND WITHOUT SURFACE SLIP CONDITIONS 
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FIG 22 VARIATION OF THE MAXIMUM TEMPERATURE WITH GAMMA 
FOR DIFFERENT VALUES OF PRANDTL NUMBER USING SUP 
AND NOSLIP SURFACE CONDITIONS 
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FIG. 23 STAGNATION POINT TEMPERATURE PROFILES WITH AND WITHOUT SURFACE 
SLIP CONDITION FOR DIFFERENT VALUES OF GAMMA 




. 24 STAGNATION POINT TANGENTIAL VELOCITY PROFILES FOR DIFFERENT 
VALUES OF GAMMA WITH SURFACE SUP CONDITIONS 










max 



Y 


FIG. 27 VARIATION OF THE MAXIMUM TEMPERATURE WITH GAMMA 

FOR DIFFERENT VALUES OF PRANDTL NUMBERS FOR SURFACE 
SLIP CONDITIONS 
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FIG. 30 

FIG. 29 VARIATION OF STAGNATION POINT SKIN-FRICTION COEFFICIENT PARAMETER, 
Zpl sin 0 , WITH GAMMA FOR NOSLIP CASE 
FIG. 30 VARIATION OF STAGNATION POINT SKIN-FRICTION COEFFICIENT PARAMETER, 
Cf/ sin 6, AGAINST GAMMA WITH AND WITHOUT SURFACE SLIP CONDITION 
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FIG. 31 COMPARISON OF STAGNATION POINT HEAT TRANSFER COEFEICIENT AT DIFFERENT VALUES 
OF GAMMA FROM PRESENT RESULTS WITH DSMC RESULTS (Ref. 9) 
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FIG. 32 COMPARISON BETWEEN STAGNATION POINT SKIN-FRICTION COEFFICIENT PARAMETER 

Cf /sine FOR DIFFERENT VALUES OF GAMMA FOR VARIOUS SPACE SHUTTLE 
RF FNTRY rONniTIONS 









FIG. 37 STAGNATION POINT TEMPERATURE AND VELOCITY PROFILES WITH AND WITHOUT 
ABLATION WITH NOSLIP SURFACE CONDITIONS 
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FIG. 39 VARIATION OF HEAT -TRANSFER COEFFICIENT, C H ,W!TH 

Mass transfer parameter ,m ,for surface slip 

AND NOSLIP CONDITIONS 
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FIG- 40 VARIATION OF THE SKIN- FRICTION COEFFICIENT PARAMETER 
Cf / sin 9 , WITH MASS TRANSFER PARAMETER M,FOR 
SURFACE SUP AND NO SLIP CONDITIONS 




